Reversal of Thermal Rectification in Quantum Systems 
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We study thermal transport in anisotropic Heisenberg spin chains using the quantum master equa- 
tion. It is found that thermal rectification changes sign when the external homogeneous magnetic 
field is varied. This reversal also occurs when the magnetic field becomes inhomogeneous. Moreover, 
we can tune the reversal of rectification by temperatures of the heat baths, the anisotropy and size 
of the spin chains. 
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Considerable progress has been made both theoreti- 
cally and experimentally on thermal transport in micro 
and nano scale in last decade [l[. It has been found that 
similar to electrons, the heat due to phonons can be 
used to carry and process information In particu- 
lar, several conceptual thermal devices have been pro- 
osed, such as thermal rectifiers Q, thermal transistors 
thermal logical gates [j|, thermal memory 0, some 
molecular level thermal machines 0, H|, and thermal 
ratchet Q. Much work has also been done in quantum 
heat transport of nanostructures an d sp in systems 
[U El El El E! EI E3, El, EI IH EH, [S| where the 
magnetic field is another degree of freedom to control 
heat flow. Indeed, it is demonstrated that thermal rec- 
tification and negative differential thermal resistance are 
observable in quantum spin chain by applying a nonuni- 
form magnetic field [231 ]. 

In this paper, we would like to concentrate on the ther- 
mal rectification in a quantum spin model. Our primary 
interest is to understand how far we can control the heat 
flow by tuning the system parameters such as magnetic 
field, system size, configurations of the system etc. In 
particular, we would like to see whether the reversal of 
thermal rectification, which has been observed in classi- 
cal system [24|. can happen in such a quantum system. 

We consider a Heisenberg spin-1 /2 chain, whose Hamil- 
tonian reads 



N-l 



N 



(i) 

where N is the number of spins, the operators of, of and 
of are the Pauli matrices for the ith spin, J x , J y and J z 
are the coupling constants between the nearest-neighbor 
spins, and hi is the magnetic field strength at the ith site. 
We set J x = J (I + 7), J y = J {I — 7), J z = J to consider 
the anisotropy in x — y plane, where 7 is the anisotropy 
parameter, and J — I, without loss of generality. Figure 
1 shows a schematic representation of this model. 
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FIG. 1: (color online) A schematic representation of the 
model with size N = 8. The spin chain is connected to two 
phonon baths with coupling Xl and Xr. The phonon baths 
are at different temperature Tl and Tr. A magnetic field is 
applied to the spin chain. 



The total Hamiltonian including two baths is 



H, 



H + H f 



Hr 



(2) 



Here Hb is the Hamiltonian of the heat baths Hb = 

fff, Hg = Ujbfbj, where 6+ and bj are 

k=l,r ' jeK 

phonon creation and annihilation operators with the 
mode cjj-; Hi is the interaction between the spin chain 

Yk, where 



3 ' 

and phonon heat baths, Hj = ^2 Xk 

K=L,R 

X L = af,Xfj = a'fj, and Yk is bath operator Yk = 



jeK 

We use the quantum master equation (QME) method 
(Refs.fH, [H, [26|, [27| El) to study heat conduction in 
this model. By tracing out the baths within the Born- 
Markovian approximation, we obtain the equation of mo- 
tion for the reduced density matrix of the system (h = 1), 



dt P 



-i[H,p] + C L p + C R p, 



(3) 



where ClP and Crp are the dissipative terms due to the 
coupling with the left and right heat bath. ClP is given 
by ClP — [XlP, Xl] + h.c. and Crp can be given in the 
similar way. Here the operator Xl can be written as 



(m\X L \n) = \L£ m .nN L {e m ,n){m\X L \n), 



(4) 



where e m ^ = e m - e n and N L {e m , n ) = (e £m -"/ TL - 
is the Bose distribution (ks = 1) with Tl being the tem- 
perature of left heat bath. \n) and e n are the eigenstates 
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FIG. 2: (color online) Rectification as a function of the mag- 
netic field in a chain of size N — 6. The temperature of the 
heat bath are T+ = T (l + A) and T_ = T (l - A), where 
To — 0.1 is the mean temperature. The coupling between 
the spin chain and heat bath are Al = 0.16 and A_r = 0.04, 
anisotropy parameter 7 — 0.1. The square, circle and trian- 
gular correspond to A = 0.2, 0.4 and 0.8, respectively. 



and eigenvalues of the spin chain system. The bath spec- 
trum function we used is of an Ohmic type. Assuming 
that the temperature is high enough to make dephasing 
fast 27], we can solve the resulting kinetic equations of 
the state probabilities numerically. The evolution time 



is chosen long enough such that the final density matrix 
reaches a steady state p s t, that is, p nn = 0. In the steady 
state, J2 n £ nf>nn = Jl + Jr = 0, then we can get the 
heat current as J = Jl = — Jr, if Tl > Tr. To quantify 
the rectification efficiency, we define rectification R as: 
R = (J + — J_)/max{ J + , J-}, where the forward heat 
flux J + is the heat current (from left to right) when the 
bath at higher temperature T + is connected to the left 
end of the chain and the backward flux J_ is the heat 
current (from right to left) when the left end of the chain 
is in contact with the bath at lower temperature T_ . 

The quantum spin chain is a nonlinear system; if we 
introduce asymmetry to the system, then it may show 
the rectification effect. We connect the spin chain to 
the two heat baths by different couplings A^ = 0.16 and 
An = 0.04. As is shown in Fig.[2l when the applied mag- 
netic field increases, the rectification changes sign. This 
phenomenon is not observed in electronic counterpart. 
From Fig. [2j we can see that rectification R can be pos- 
itive, zero, or negative, depending on the magnetic field 
h. The behavior of R remains similar for different A, the 
temperature difference of the baths. 

The mechanism of thermal rectification reversal can be 
understood from a two-spin system which has four eigen- 
states, whose contribution to heat transport can be seen 
clearly in Fig. [3l Fig. [3ja) shows the eigenvalues E n vs 
magnetic field h, where two eigenvalues do not change 
with the magnetic field. Fig. E[b) shows the steady den- 
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FIG. 3: (color online) Analysis of the reversal of rectifica- 
tion in a two-spin system. Parameters: To = 1.0, A = 0.8, 
\l = 0.16, \r — 0.04, 7 = 0.1. (a) The eigenvalues E n vs 
magnetic field h. (b) The density p nn (probability of each 
eigenstate) vs h. (c) The product of E n and the derivative of 
p nn from high temperature heat bath, which is the contribu- 
tion to heat flux from each eigenstate. In (a), (b), and (c), 
the square, circle, up triangle, and down triangle correspond 
to eigenstate 1, 2, 3, and 4, respectively. In (b) and (c), the 
solid and hollow symbols correspond to forward and backward 
thermal transport, respectively, (d) The heat fluxes J+, J_ 
vs magnetic field h (left scale). The square curve shows rec- 
tification effect R vs magnetic field h (right scale). 



sity p nn as a function of magnetic field. From this figure, 
we find that when the magnetic field is weak, the ground 
state and the excited states have some probability to be 
occupied (the ground state has the largest probability) . If 
the magnetic field increases, the probability of the ground 
state becomes larger; after a certain value, the probabil- 
ity of ground state is close to one, while others are zero. 
The heat current from the contribution of each eigenstate 
is shown in Fig. [3fc) . When the magnetic field is weak, 
each state contributes to the heat flux either positively 
or negatively. We find that the contribution from ground 
state in the backward thermal transport is larger than 
the forward case. But the contributions from other ex- 
cited states have negative effect; therefore the total effect 
is the forward heat flux larger than the backward flux, 
that is, the rectification is positive. When the magnetic 
field increases over a certain value, the ground state will 
dominate the heat transport, and the contribution from 
other states decrease. At the time the total heat flux will 
have the similar behavior with the contribution of ground 
state: the backward flux is larger than the forward one; 
both of them increase first and decrease to zero at last, 
when the system will stay in the ground state. There- 
fore, the rectification changes from positive to negative, 
then from negative to zero at last, which can be seen in 
Fig- Eld). In short, the heat flux from the contribution 
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FIG. 4: (color online) The rectification R vs nonuniform mag- 
netic field. The chain size is TV = 6, h(l : N/2) is the mag- 
netic field applied to site 1 to TV/2, and others are zero, that is, 
h{N/2 + l : TV) = 0. Here 7 = 0, T = 0.1 and \l = An = 0.1. 
The square, circle and triangular correspond to A = 0.2, 0.4 
and 0.8 respectively. 
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FIG. 5: (color online) Rectification as a function of mean tem- 
perature To. The size of spin chain is N = 6. The temperature 
of the heat bath are T+ = To + AT and T_ = To - AT. Here, 
\ L = 0.16, \ R = 0.04, 7 = 0.1 and h = 0.01. The square, 
circle and triangular correspond to AT = 0.01, 0.02 and 0.04 
respectively. 



of ground state is larger when the more weakly coupled 
reservoir is hotter, that is, the backward one is larger 
than the forward one; however for the flux from excited 
states, the forward one is larger than the backward one. 
The rectification is determined by the competition of the 
contributions to heat flux from the ground state and the 
excited states. 

Applying a nonuniform magnetic field to the spin chain 
is another possible way to introduce asymmetry to the 
system, and the system can also exhibit rectification. 
Fig. [4] demonstrates this phenomenon in such a system. 
The rectification is zero when the magnetic field is zero 
because of no asymmetry in the system. When a weak 
magnetic field is applied to the left half part of the spin 
chain, the rectification becomes positive. When mag- 
netic field increases, the energy difference increases for 
the left part, which enlarges the rectification. If the field 
increases further, then the rectification reverses. For dif- 
ferent temperature difference, the rectification has simi- 
lar effect but different magnitude. 

From the above discussions, we find that the rectifi- 
cation can change sign with the applied magnetic field. 
Indeed, an increase of mean temperature can also induce 
reversal of rectification. Fig. [5] shows that the rectifica- 
tion as a function of mean temperature To. Here we keep 
the temperature difference fixed, and increase the mean 
temperature of the heat baths; the rectification changes 
sign from negative to positive. When the temperature is 
very low, the ground state dominates the thermal trans- 
port, the backward flux is larger than the forward one; 
if we raise the temperature, more excited states will con- 
tribute to the heat flux and gradually control the thermal 
transport, when the rectification changes to positive. For 
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FIG. 6: (color online) The rectification changes with magnetic 
field for the spin chain with different sizes. The curves of 
square, circle and up triangle, and down triangle, star and 
diamond correspond to N = 2, 4, 6, 8, 10 and 12 respectively. 
Here the parameters are: T+ = 0.18 and T_ = 0.02, Xl — 
0.16, \r = 0.04 and 7 = 0.1. 



different temperature difference, the rectification changes 
sign at almost the same mean temperature. 

Rectification can change sign with the external param- 
eter, such as magnetic field and the temperatures of the 
heat baths. In our study, we find that thermal rectifica- 
tion can reverse with the properties of the spin chain it- 
self, such as the size and the anisotropy of the spin chain. 
Fig. [6] shows the rectification changes with the magnetic 
field for different size cases. The rectification effect be- 
haves differently for different size. In Fig. there is no 
reversal of rectification for small size cases TV = 2 and 
N = 4; but for larger size cases N = 6, 8, 10 and 12, it 
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FIG. 7: (color online) The rectification changes with 
anisotropy in the spin chain. The chain size is N = 6. The 
square, circle and triangular correspond to A = 0.2, 0.4 and 
0.8 respectively. Here the parameters are: T+ = 0.18 and 
T- = 0.02, \ L = 0.16, Xr = 0.04 and h = 0.01. 



shows reversal of rectification. Fig. [7] shows the rectifi- 
cation reverses with the anisotropy of the spin chain. In 
the weak anisotropy range, the rectification coefficient is 
positive, but it changes to negative when the anisotropy 
is strong. During the changing of anisotropy, the for- 
ward total flux is larger than the backward one at first, 
and then reverses, although the heat contribution from 
ground state in the backward transport is always larger 
than that in the forward one. 

In conclusion, we have studied thermal rectification 
in quantum spin-chain systems by using quantum mas- 
ter equations. It has been shown that rectification 
can change sign when the magnetic field, temperature, 
the anisotropy, and the system size change. Although 
the reversal of rectification is complicated parameter- 
dependent, it is believed to be a universal phenomenon 
for the thermal transport in one-dimension systems. 
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